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INTRODUCTION 

There  have  been  a  number  of  theoretical  studies  of  the  amplification 


of  sound  propagating  through  eases  under  a  variety  of  non-eou i 1 i hr i  uni 


1-3 


conditions.  Tngard  and  others  have  nredicted  the  amplification  of 

i-h 


sound  nronayating  in  a  weaklv  ionized  gas.  Cilhert  and  others 


predicted  the  ampl  i  f  i  ca  t  ion  of  sound  propagating  in  reactin'1 


7  8  q  in 

Plomb,  Pchulz,  Rauer  and  Hass,  and  Shields  predicted  the  amplifica¬ 


tion  of  sound  propagating  in  gases  subjected  to  externa]  energy  input 
(energy  pumping-).  In  the  latter  cases  two  different  mechanisms  leading 
to  the  amplification  of  sound  have  been  proposed:  selective  absorption  of 
the  oumping  energy  and  relaxation  processes  of  internal  modes.  The  sound 
amplification  due  to  the  relaxation  processes  results  from  tKe  Faet  that 
the  vibrational  energy  transition  rate  is  pressure  and.  temperature  depen¬ 
dent.  The  increase  in  pressure  and  temperature  results  in  an  increase  in 
the  transition  rate,  which  in  turn  results  in  a  further  increase  in 
temperature  (here,  temperature  means  translational  temperature).  THns, 
in  the  presence  of  the  sound  wave,  the  vibrational  relaxing  energy  will 
he  selectively  transferred  into  the  hot  part  of  the  sound  cycles.  Exper¬ 


imental  studv  to  verify  these  theories  is  still  limited.  Recently  Hass 


1  1 


and  Detsch  observe'!  the  amplification  hv  gas  phase  reactions  in  Cl-**  - 


\r  and  Cl-1'  -SF  mixtures.  In  this  studv,  the  amplification  of  sound  hv 

>  h 


vibrational  relaxation 


s  is  of  primary  interest. 


T  p  p  pH  p  f  to 


1 


? 


separate  the  effect  due  to  relaxation  from  that  due  to  pumping,  the 
pumping  mechanism  is  terminated  before  the  sound  wave  is  introduced  into 
the  gas. 

When  energy  is  pumped  into  the  gas  and  then  the  numoing  process 
is  terminated ,  a  metastable  state  with  an  overpopulation  of  excited  vibra¬ 
tional  states  results.  Relaxation  to  equilibrium  occurs  due  to 
vibratonal-translational  (v-t)  energy  transfer,  or,  in  the  case  of  polya¬ 
tomic  gases  due  to  vibrational-vibrational (v-v) ,  followed  bv  the  v-t 
energy  transfer.  The  amplification  of  sound  wave  should  occur  when  it 
is  propagated  through  the  gas  during  the  relaxation  process.  In  this 
case,  however,  the  sound  period  must  be  short  compared  to  the  relaxation 
time.  Shields^  has  considered  three  typical  situations  in  which  it  is 
expeimentally  feasible  to  observe  the  amplification  of  sound.  The  three 
situations  are: 

1.  Metastable  state  with  a  single  long  lived  relaxing  mode: 

2.  Metastable  state  with  fast  v-v,  slow  v-t  exchanges;  and 

3.  Metastable  state  with  slow  v-v,  fast  v-t  exchanges. 

The  first  case  can  be  produced  in  most  diatomic  gases.  Diatomic  gases 
such  as  N^,  C>2 ,  Cl^,  and  CO  have  a  single  vibrational  mode  with  a  long 
relaxation  time.  After  energy  is  pumped  into  excited  vibrational  levels 
a  metastable  state  persists  for  a  period  of  time  longer  than  the  period 
of  low-f renuencv  sound.  The  last  two  cases  can  be  produced  in  polvatomic 
gases  or  some  mixtures.  Polvatomic  molecules  possess  several  vibrational 


3 


modes  which,  in  principle,  can  relax  independently.  Most  polyatomic 

12  13 

gases,  however,  show  a  single  relaxation  time.  ’  This  is  due  to  the 

rapid  v-v  energy  transfer  which  maintains  continuous  eauilibrium  between 

modes.  The  amplification  of  sound  propagating  in  such  gases  could  be 

observed  if  the  sound  period  is  short  compared  to  the  v-t  relaxation 

time.  In  some  gases  the  v-v  transfer  may  be  slower  than  the  v-t 
12  14 

transfer.  ’  In  this  situation  the  sound  period  must  be  short  compared 
to  the  v-v  relaxation  time.  As  examples  of  the  last  two  cases,  Shields1 2^ 
calculated  the  expected  sound  amplification  in  CD ^  and  in  a  N9-He-CO^ 
mixture. 

In  this  investigation  the  sound  wave  is  propagated  through  a  vibra- 
tionallv  excited  pure  nitrogen  gas  and  mixtures  of  nitrogen  with  helium 
and  hvdrogen.  The  overpopulation  of  vibrationallv  excited  states  is 
produced  by  an  electrical  discharge.  A  sound  pulse  is  generated  in  the 
gas  by  the  sudden  increase  in  gas  temperature  accomnanving  the  discharge. 
In  this  situation  the  sound  pulse  propagates  immediately  after  the  meta¬ 
stable  state  has  been  established,  so  that  the  requirement  discussed 
above  is  satisfied.  Nitrogen  gas  was  chosen  for  the  following  reasons: 

1.  Nitrogen  molecules  have  long  relaxation  times  of  the  order  of 
seconds  at  atmospheric  pressure;  and 

2.  Studies  have  shown  that  a  large  percentage  of  the  energy  from 
the  electrical  discharge  will  wind  up  in  vibrational  energy  of 


nitrogen  molecules 


15 


\'itrogen-hel  ium  and  ni C rogen-h vd rogen  mixtures  were  used  to  see  the 

effects  of  changing  the  vibrational  relaxation  time  and  the  thermal 

conduction  relaxation  time.  The  addition  of  helium  and  hvdrotren  mole- 

1  h-1  Q 

cules  drastically  decreases  the  relaxation  time  of  nitrogen. 

The  sound  absorntion  and  velocity  in  the  gas  were  determined  bv 
observing  the  decav  and  the  transit  time  of  the  comnressional  Dulse  that 
reflected  back  and  forth  in  the  gas.  To  determine  the  variation  of  the 
absorption  and  velocity  (the  translational  temperature  was  determined 
from  the  velocity)  with  frequency  and  time,  the  nulse  was  Fourier 
analyzed  after  traveling  various  distances  in  the  gas  sample,  and  the 
absorption  and  velocity  were  determined  for  each  frequency  component. 


CH'VPTFP  TT 


theory 


CALCULATION!  OF  THE  SOHNn  ABSORPTION  AMT)  VELOCITY  IN  AM  POUTLT P.P^'M  CAS 

As  discussed  in  the  introduction,  the  amplification  of  sound  bv 
relaxation  nrocesses  results  from  the  preferential  heating  of  the 
compression  part  of  the  sound  wave  from  the  vibrational  relaxation 
energy.  However,  the  gain  in  energy  will  he  countered  Hv  losses  due  to  a 
number  of  effects.  To  observe  the  effects  of  the  excited  vibrational 
modes  on  the  sound  propagation  it  is  therefore  necessarv  to  determine  the 
absorption  and  velocity  in  the  unexcited  gases. 

The  most  important  effects  leading  to  the  sound  absorption  are  shear 
viscositv,  heat  conduction,  loss  to  the  wall  of  the  tube,  and  thermal 
relaxation.  The  absorption  resulting  from  viscosity  and  thermal  conduc¬ 
tion  is  termed  classical  absorption.  The  classical  eauation  for  sound 

oi 

absorption  and  velocity  in  a  tube  was  first  developed  bv  Firchoff. 

Beginning  with  the  conservation  eouations  for  mass,  energy,  and  momentum, 

he  obtained  an  algebraic  enuation  for  the  propagation  constant  for 

radially  symmetric  waves  bv  approximating  Bessel  functions  with  a  one 

term  expansion.  In  this  derivation  he  assumed  that  the  particle  velocity 

22 

and  the  sound  temperature  were  zero  at  the  tube  wall.  Shields  Has 
solved  for  the  absorption  and  velocity  numerically  and  included  the  slin 
velocity  and  the  temperature  "jump  at  the  tube  wall.  In  this  study  the 


.%  \  .S  V 
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sound  absorption  and  velocitv  in  the  un excited  gas  were  calculated  at 

various  freauencies  and  temperatures  using  the  numerical  solution 

23 

developed  bv  Shields.  The  computer  program  for  this  numerical  solution 

allows  for  calculation  of  viscous  and  thermal  losses  both  in  the  bodv  of 

the  gas  and  at  the  tube  wall  and  for  relaxation  effects.  For  this 

purnose  the  viscositv  of  erases  used  in  this  studv  were  computed  over  the 

temperature  ran?e  from  290 aK  to  350°Y  using  formulae  and  potential 

parameters  of  the  gas  given  in  Reference  23.  The  specific  heat  for  pure 

24 

gases  were  taken  from  the  Table  of  Thermal  Properties  of  Oases.  The 

specific  heat  of  mixtures  were  obtained  bv  taking  a  weighted  average  of 
the  pure  component  specific  heats.  The  relaxation  times  were  taken  from 
the  work  of  Frey,  et  al.  ^  and  Audibert,  et  al.'^.  the  momentum  accommo¬ 
dation  coefficient  and  the  energy  accommodation  coefficient  were  set 
eaual  to  1.0  and  0.8  respectively.  The  values  are  approximate  but  have 
little  effect  on  the  calculated  absorption  and  velocitv.  The  physical 
properties  of  gases  under  studv  are  given  in  Appendix  R.  To  check  the 
accuracy  of  this  calculation,  measurements  on  stable  pure  nitrogen  gas 
and  helium-nitrogen  mixtures  were  made  at  various  pressures  and  freauen¬ 
cies.  The  experimental  method  and  results  are  discussed  in  Chanters  ITT 
and  TV. 

CALCULATION  OF  THE  ABSORPTION  AND  VELOCITY  IN  TV>E  CAS  UT T”  A  mfta^TARLE 
STATE 

In  this  section  the  eauations  leading  to  the  amplification  of  sound 
are  developed,  and  the  numerical  solutions  to  the  resulting  enuations  for 


several  gases  are  presented 


Assuming  that  the  attenuation  and  gain 


rates  are  additive,  the  amplification  factor  mav  he  determined 
separatelv. 

The  relaxation  absorption  (or  amplification  as  the  case  mav  he)  and 
velocitv  can  be  obtained  from  the  propagation  constant  for  the  acoustic 
wave.  Using  the  continuitv  equations  for  mass,  momentum  and  energv  and 
neglecting  the  internal  friction  and  thermal  conduction,  the  propagation 
constant  is  given  bv 

k  =  /V/RT  /3p/3P.  (1) 

For  an  ideal  gas  Fauation  1  becomes 

k  =,Vm/rt  /c  /C  ,  (2) 

V  p 

where  ui  is  the  angular  frequencv,  M  is  the  molecular  weight,  P  is  the  mas 

constant,  T  is  the  absolute  temperature,  and  C  and  C  are  the  specific 

v  p 

heats  at  constant  volume  and  pressure.  For  a  relaxaing  gas  C  and  C  are 

d  v 

comolex  and  frequencv  dependent.  Writing  k  in  Equation  2  in  the  complex 
form 


k  =  kj  -  ik2, 


(2) 


2 

where  i  =  -1 ,  the  sound  absorption  (a  )  and  the  velocitv  (v  )  are  given 

s  s 

by 


1  W 

ct  =  k  ;  v  =  : — . 
s  2  s  '  k 


(4) 


1- 


Sound  amplification  results  if  is  negative.  The  problem  of 

calculating  the  absorption  and  velocity  is  now  reduced  to  finding  the 

real  and  imaginary  part  of  C  and  0  which  can  he  derived  from  the 

p  v 

reaction  eauation  that  gives  the  time  rate  of  change  of  energy  in  a 

in 

vibrational  mode.  The  following  was  adapted  from  the  work  of  Shields. 

When  nitrogen  gas  is  subiected  to  an  electrical  discharge  most  of 

the  energy  deposited  in  the  gas  will  wind  up  in  the  vibrational  energy  of 

molecules, ^  leaving  the  gas  in  the  state  with  an  overpopulation  of 

excited  vibrational  states.  Tiue  to  molecular  collisions  some  of  the 

vibrational  energy  will  be  converted  back  to  the  translational  energy  bv 

v-v  and/or  v-t  exchanges.  A  sound  wave  propagated  during  this  transition 

process  will  cause  small  sinusoidal  perturbations  in  the  slowlv  decaving 

values  of  temperature  and  pressure.  Assuming  that  the  normal  modes  of 

molecular  vibration  can  be  treated  as  harmonic  oscillators,  and  that 

resonant  exchange  of  energy  between  levels  within  a  single  mode  rapidlv 

establishes  a  Boltzman  distribution,  the  rate  of  change  of  energy  in  the 

25 

ath  vibrational  mode  can  be  written  as 

.  a^-a+i ;  8->-g  +  i  .  .  na(a)n^(b) 

Ea  ■  ih  lhV,£oeSo  ,k  u’b)  —5; - 

Nb 

where  a  and  g  refer  to  the  vibrational  Quantum  number  of  modes  a  and  b 

respectively,  P  is  the  pressure,  N  is  the  total  number  of  molecules 

b 

carrying  mode  b,  n  and  n  are  the  numbers  of  molecules  occupying  and 


levels  in  modes  a  and  b  respectively,  h  is  the  Planck's  constant,  is 


carrving  node  h,  na  and  ng  are  the  numbers  of  molecules  occunvin?  a  and  3 


r 

K 


levels  In  nodes  a  and  b  respectively,  h  is  the  Planck’s  constant,  v  is 

a 

a-*a+i 

the  fundamental  vibrational  frequency  of  mode  a,  and  k  is  the 

composition  average  rate  constant  in  the  dimension  (time  pressure)  ^  for 

the  particular  process  in  which  a  changes  to  a+i  and  3  to  3+i.  The 

summation  over  a  and  3  gives  a  net  rate  constant  for  transitions  where 

mode  a  gains  i  quanta  and  mode  h  gains  i  auanta. 

Equation  5  assumes  that  mode  a  is  exchanging  energy  with  mode  b  and 

with  translation  ( i  =0 ,  n  /U  =1).  Modes  a  and  b  can  be  in  the  same  mole- 

b 

cules,  or,  in  the  case  of  mixture,  in  different  molecules  and  mav  have 
different  temperatures.  Assumptions  that  the  normal  vibrations  can  be 
treated  as  harmonics  oscillators  and  a  Roltzman  distribution  is  main¬ 
tained  allows  us  to  write  the  following  relations 


ka^a+i,3+6  +  i  „  Ha(a-l)(a-2) _ (a-i+Dl, 


x  f  (B+i )  (B+i-1 ) . .  • .  (3  +  0  1  /i  !  j  M  x  k 


i-*-n ,  o-»-i 


(f 


and 


/A  t  \  ~  hva/kTa  , ,  -hva/k._.  -  hv_/kTa  ,  7 . 

nr. (a)  =  n.(a)e  a  a  =  N  ( 1  — e  a  a)e  a  a,  (/) 

u  0  a 


where  M  is  the  number  of  molecules  with  node  a.  Using  Equations  <S  and  7 
a 


and  the  relation 


(ct+i)(a  +  i-l)....(a4i)x'-‘  =  -^(7^—) 


X! 


dx1  1-x  (1+x) 


i+1 


(3) 


Equation  5  can  be  written  as 


1° 


E  =  .  Z  .  i  N  hv  P(kr"0,0'M  (a,b)e"10a/Ta 
a  i  ,  i  a  a 

_  kOi  ,i^O(a>b)e-i0b/Tb)/(fx  (T  )]ifx  (T 


a  a 


kb  h 


)i). 


(9) 


hv 


^  — 0  j  'j* 

where  0  =  — r — ;  y  (T  )  =  1-e  a  a.  '’ere  the  principle  of  detailed 

a  *  a  a 

i-H+i  j  +i-M  -iQ /T 

balancing,  k  =  k  e  ,  has  been  used.  Enuation  9  can  now  be 

specialized  to  the  three  situations  discussed  in  the  introduction. 

Most  diatomic  eases  have  onlv  a  single  vibrational  mode  with  a  lone 
relaxation  time.  In  this  situation,  the  vibrational  mode,  sav  mode  a, 
communicates  onlv  wi*-h  the  translational  mode  so  j=0  and  i=l.  Equation  9 
then  becomes,  after  some  manipulations, 


E  =  -Pk10v  (T  )  f  F.  (T  )-F.  ( T>  ] 
a  Aa  a  a  a  a 


(10) 


.  ,  Nak9aexP(-0a/Ta)  ,  ,  Nak9aexp(-0a/T) 

“here  W  *  - fl7la~  a"d  E»‘T)  ‘  xa(l)  “he"  a 

sound  wave  is  introduced  F  will  depart  from  its  steadv  state  value.  Let 

a 

•  •  • 

AF.  be  the  change  in  E  due  to  the  presence  of  the  sound  wave.  AF  can 
a  a  a 

be  obtained  bv  differentiating  Eauation  9  with  resnect  to  the  transla¬ 
tional  temnrature  T.  The  result  is  given  bv 


,  r  1  dP  ,  1  dk10  9aexp(ea/T)  AEa(Ta)-AEa(T)  ,ttx 

=  F-  H  p  TT  +  TrT  Du? - 5 - 1 - 1AT - - - )  *  U  1  ) 

a  a  P  dT  k10  dT  Dxa(T) 


^2 


Here  onlv  the  first  order  terns  of  the  sound  variable  are  maintained. 


f1-,  T  ,  and  T  represent  values  without  the  sound  wave.  Introducing  the 
a  a 

sinusoidal  time  dependence  and  setting  AF  (T  )  =  C  (T  )AT  ,  AF(T)  = 

a  a  a  a  a 

,  .  ‘0  Ea(Ta)-Ea(T) 

C  (DAT  and  F  =  -  — - — -  enuation  11  becomes 

a  a  t 


11 


ATa  E§t  i  hp 

(l+iWT)“=  [Ca(Ta)]  X  [P  dT  + 

Ca(T) 

'  Ca(Ta) 


d ( Ink 10 (a) ) 


eae-9*/T 

T2Xa(T) 


In  the  case  in  which  the  gas  is  in  equilibrium  (the  gas  has  not  been 

pumped)  T  =T  and  E°=ft.  Equation  12  becomes 
a  a 

AT,  -1 

“  =  ( l+iu>T  )  (13) 


and  the  frequency  dependent  specific  heat  is  given  bv 


AT  a 

C>>  = 


Eoo  +  7—1 -  ;  C  (ai)  =  C  (oj)+R 

1+j  WT  O  v 


where  C®  is  the  frequency  independent  part  of  the  specific  heat.  When 

the  gas  has  been  pumped  the  vibrational  temperature  is  not  eaual  to 

the  translational  temperature  T  and  E"  is  not  zero.  The  frequency  depen- 

a 

dent  specific  heat  is  now  given  bv 

,  E°x 

C  (w)  =  C*,  +  C’O+jwT)  x;  C  (w)  =  C  (w)+R  -  z— - - r,  (15) 

v  a  p  v  T(1+jcjt) 

where 


Ca  =  [1  +[d-lnk--(--)]--^ - ]  +  C  (T  )  (16) 

a  T  dT  Ta2xa(T)  3  3  , 

and R is  the  universal  gas  constant.  The  extra  term,  E°/T(l+jwT),  the 

a 

second  equation  of  Equation  14  results  from  the  fact  that,  at  constant 
9P 

pressure,  -n;  =0.  Amplification  will  occur  if  C'  becomes  negative.  This 
3T  a 

is  the  case  when  T  -T  is  a  few  tens  of  degrees  and  d(lnk^)/dT  (with 

a 

dimension  K  ')  is  of  the  order  of  a  few  one  hundreds.  however,  the 


amplification  may  be  small  and  too  difficult  to  observe. 


Hsing  Equations 


2,4,15,  and  16,  Shields  has  calculated  the  amnli f ication  factor 


several  diatomic  gases.  The  results  are  shown  in  Table  1. 


TABLE  1 

Sound  amplification  in  dB  per  wavelength  for  tvpical  diatomic  gases 
at  300°K.  d(lnK  )/dT  is  set  equal  to  0.01,  and  ml  =100. 


GAS 

Amplification 

factor 

(dB/A) 

(cm)  1 

T  =500°K 
a 

|  1000  “K 

1 

2000  °K 

|  3000° 
1 

K  |  4000  °K 

1 

1  6onn°K 
! 

N2 

233  1 

0.0027 

(  0.082 

0.52 

1  1.1 
| 

|  1.7 

| 

|  2.9 

i 

CO 

2143 

0.0042 

1  n. 10 

1 

0.56 

1  1.2 

1 

|  1.8 
| 

|  3.0 

1 

°2 

1556 

0.015 

1 

|  0.18 

0.73 

1  1.4 

i 

1 

1  2.0 

1 

1 

I 

557 

0.075 

|  0.37 

1.00 

1 

! 

1 

i 

i 

In  this  calculation  ojt  was  set  eaual  to  100  and  -^-(lnk'0)  equal  to  0.01. 
For  this  range  of  ujT  the  amplification  per  wavelength  is  approxi matelv 
proportional  to  (wx )  ^  or  the  amplification  per  unit  length  is  indepen¬ 
dent  of  frequencv. 

The  equations  derived  above  also  apolv  to  a  mixture  of  a  diatomic 
gas  with  a  monoatomic  gas  such  as  N^-He.  In  pure  the  relaxation  time 
is  of  the  order  of  atmosphere  seconds.  To  observe  appreciable  gain, 
measurements  would  have  to  be  made  at  such  low  frequencies  and  long 
wavelengths  that  it  would  be  verv  difficult  to  conduct  the  needed 
measurements  in  the  laboratorv.  Bv  adding  He  or  H0  molecules,  t^e  Ire- 


13 


When  the  sound  wave  is  propagated  through  a  metastable  state  with 
more  than  one  excited  vibrational  mode  the  oroblem  becomes  more  compli¬ 
cated.  However,  the  change  in  E^  due  to  the  sound  wave  perturbation  can 
be  obtained  bv  differentiating  Equation  5.  The  freauencv  dependent 
specific  heat  similar  to  Equation  14  is  given  by 


C  (w) 


Coo  +  C 

a 


£a 

AT 


+  C, 


ATb 

AT 


+ 


•  •  •  •  y 


(17) 


where 


(IbVVW 

C.  =  R  — - - - - -y 

i=a,b, .  Using  Eauation  17  and  a  similar 

has  ca1culated  the  sound  amplification  in  CD, 

4 

The  results  are  shown  in  Table  2  and  Figure  1 


(1*) 

eauation  for  C  (w),  Shields 
n 

and  a  mixture  of  N^-CO^-He. 
of  Reference  10. 


EXPERIMENTAL  CONDITION  FOR  OBSERVING  THE  AMPLIFICATION 

As  discussed  in  the  preceding  section,  an  electrical  discharge  in 

low  pressure  N  will  result  in  overpopulation  of  excited  vibrational 

states.  The  electric  circuit  used  in  this  experiment  to  produce  the 

discharge  is  able  to  deposit  approximate lv  2.5  ioules  of  energy  in  the 

28 

gas.  Assuming  about  70%  of  this  goes  to  vibrational  energy,^  the 
vibrational  temperature  was  estimated  to  he  of  the  order  of  a  few  thou¬ 
sand  degrees  Kelvin  depending  on  the  volume  and  pressure  of  the  eas. 

After  the  discharge  has  been  turned  off,  according  to  Margottin 

28 

Waclou,  Dovennette,  and  Henrv  ,  the  vibrat i onal 1 v  excited  gas  will  relax 
in  three  different  wavs.  ^irst  is  relaxation  due  to  collisional 


de-excitation  in  which  the  vibrational  energv  is  transferred  into  trans¬ 


lational  energv.  Second  is  relaxation  due  to  wall  collision  in  which  the 

vibrational  energy  is  transferred  to  the  wall.  Third  is  relaxation  due 

to  spontaneous  radiation.  Since  N0  has  no  dipole  moment,  the  latter 

process  is  not  involved  in  its  relaxation. 

Let  x  and  T  be  the  relaxation  times  for  transferrins  vi Kra t i on  a  1 
v  vc 

energy  to  translation  and  to  the  tube  wall,  and  let  t  he  the  relaxation 

tc 

time  eauating  the  gas  translational  temperature  with  the  wall  temperature 

hv  thermal  conduction.  The  variation  of  the  translational  temperature 

26 

with  time  then  depends  upon  t  ,  T  ,  and  x  .  Margott in-Mac  1  on ,  et  al. 

V  VC  t  c 

found  that  in  carbon  monoxide  x  =t  •  We  will  assume  that  the  two  are 

vc  tc 

approximately  equal  in  also. 

For  the  experiment  proposed  here,  it  will  be  necessary  for  the  sound 

period,  P,  to  be  less  than  either  T  ,  x  ,  or  t  .  Unless  P  is  less  than 

v  vc  tc 

x  the  propagation  will  be  isothermal  and  no  relaxation  effects  will  be 
tc 

observed.  If  P  is  not  less  than  x  and  x  the  metastable  state  will  not 

v  vc 

persist  over  a  number  of  sound  vibrations.  However,  if  P/x  is  small  v.t 

v  v 

will  be  large  and  the  gain  per  wavelength  will  be  small  (Notice  the  gain 

10. 

per  wavelength  Is  inversely  proportional  to  oox  ).  In  addition,  the 

total  oath  length  for  the  sound  in  the  gas,  and  thus  the  total  gain,  is 

restricted  bv  the  condition  that  the  transit  time  for  the  sound  wave 

through  the  gas  must  he  small  compared  to  the  vibrational  relaxation 

time,  x  ,  and/or  x  .  Otherwise  the  overpopulation  of  excited 
v  vc 

vibrational  states  decavs  awav  before  the  sound  wave  can  traverse  the 
gas . 


A  final  constraint  must  be  placed  upon  the  ratio  i  -x  .  ''nless 

v  t  c 

this  ratio  is  considerably  more  than  one,  hot  spots  are  likel  v  to  develop 

in  the  »as  and  there  is  not  time  for  spatial  variations  in  temperature  to 

he  smoothed  out  bv  conduction  before  the  vibrational  excitation  dies 

a wav.  What  the  above  means  is  that  there  is  a  verv  limited  ranee  of 

values  of  x  for  which  it  will  he  possible  to  observe  the  expected 

v  tc 

gain. 


Table  2  shows  the  value  of  t  and  i  (in  milliseconds’)  for  ^ases 

v  tc 

used  in  tKis  studv.  Values  of  x  were  taken  from  the  work  of  Frev,  et 

v 

16  26  ?7 

al.  and  i  's  were  calculated  using  the  formula 
tc 


1 


( 1 0 ) 


where  X  =  — is  the  coefficient  of  diffusivitv  (cm's  ),  and  (aCl  is  the 
J  L 

V 

first  zero  of  the  Bessel  function  .J^(aO  where  a  is  the  tube  radius.  Tn 


this  calculation,  the  thermal  conductivitv  ,k,  for  pure  nitrogen  was 


taken  from  the  Table  of  Thermal  Properties  of  Cases,  National  Bureau  of 


Standards,  Circular  564,  Vov.  1955,  and  for  mixtures  were  calculated 

using  equation  8.2-40  of  Hirschfelder ,  Curtis,  and  Bird  .  t  and  T 

V  t  c 


were  calculated  at  pressure  40  Torr  in  a  0.615  cm  radius  tube 


TABLE  2 


Vibrational  relaxation  time  x  and  translational 

relatation  time  x  (both  in  milliseconds)  for 
t  c 

several  gases ,  calculated  at  pressure  of  AO  Torr 
in  a  tube  of  0.615  cm  in  radius. 


Gas 

t  (ms) 

V 

X  _  ( ms ) 
tc 

Pure 

1 9x10  ^ 

11.56 

25Fe-75N2 

520 

7.16 

50He-50N2 

260 

3.79 

10H2-90N2 

52.8 

11.1b 

20H2-80N2 

26.4 

8.78 

As  we  can  see  from  Eauation  19,  x  is  directly  proportional  to  pressure 

tc 

and  inversely  proportional  to  the  snuared  tube  radius.  Since  x  is 

v 

inversely  proportional  to  pressure,  the  desired  ratio  of  x  lx  can  be 

v  tc 

obtained  by  choosing  the  proper  tube  radius  and  gas  pressure.  However, 
there  is  a  high  pressure  limit  due  to  the  difficulty  in  getting  a  uniform 
electrical  discharge  in  a  gas  above  a  few  tens  of  Torr  pressure. 
Increasing  the  pressure  also  increases  the  energy  needed  in  the  discharge 
to  produce  the  high  vibrational  temperature. 


£ 


!• 
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CHAPTER  III 

EXPERIMENTAL 


APPARATUS 

The  apparatus  used  in  making  the  sound  absorption  and  velocitv 
measurements  is  diagrammed  in  Figure  1.  The  different  components  are 
described  below. 

Sound  Tube 

The  sound  tube,  which  contained  a  movable  reflector,  was  a  ovrex 

tube,  about  100  cm  long.  Two  different  size  tubes,  2.54  cm  and  1.23  cm 

in  inner  diameter,  were  used.  One  end  of  the  glass  tube  was  i oined 

coaxially  to  a  short  section  of  a  brass  tube.  A  receiving  microphone  was 

mounted  on  the  side  and  a  stationary  reflector  was  mounted  in  the  end  of 

the  brass  section.  (In  experiments  with  stable  gases  the  stationarv 

reflector  was  replaced  by  a  speaker  which  served  also  as  a  reflector.) 

The  brass  section  was  connected  to  the  gas  handling  system  by  means  of 

0.25  inch  in  diameter  copper  tubing.  At  the  iunction  between  the  glass 

tube  and  the  brass  section  a  metal  screen  made  of  a  piece  of  cooper  wire 

2 

mesh  of  about  25  cm~  was  mounted  terminating  the  electrical  discharge. 
The  brass  section  together  with  the  gas  handling  system  were  grounded. 
The  other  end  of  the  tube  was  terminated  with  the  second  brass  section, 
the  end  of  which  contained  an  O-ring  seal  through  which  the  rod  moving 
the  reflector  passed.  This  second  brass  section  was  connected  to  the  gas 
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FIGURE  1 


HI AGRAM  OF  MEASURING  SYSTEM 


handling  svstem  hv  means  of  a  3/8  inch  pvrex  tube,  and  the  reflector  was 
connected  to  the  nositive  terminal  of  the  discharge  circuit.  This  end  of 
the  sound  tube  was  electricallv  isolated  from  the  ^rounded  end  bv  means 
of  plastic  rod  supports.  The  movable  reflector  and  the  screen  served  as 
a  pair  of  electrodes  between  which  the  discharge  occurred.  The  gas  in 
the  region  between  the  electrodes  was  ranidlv  heated  by  the  discharge 
whereas  the  gas  in  the  region  between  the  screen  and  the  statinnarv 
reflector  remained  cool.  The  length  of  this  cool  section  was  adjusted 
producing  the  optimum  sound  pulse. 

Microphone  and  Speaker 

The  microphone  was  mounted  on  the  side  of  the  grounded  brass  sec¬ 
tion,  apnroximateiv  3.5  cm  awav  from  the  screen,  and  immediatelv  in  front 
of  the  stationarv  reflector.  Its  face  was  parallel  to  the  tube  axis  and 
was  adjusted  so  that  the  microphone  surface  protruded  onlv  slightlv  into 
the  sound  tube  avoiding  sound  reflections.  In  experiments  with  the  2.54 
cm  tube,  the  microphone  was  an  adapted  one  inch  Bruel-Kiaer  (B-K)  conden¬ 
ser  microphone.  The  microphone  was  connected  to  the  nreamp  through  a 
vacuum  feed-through  that  was  sealed  through  the  wall  of  the  brass  tube 
with  a  Swagelock  fitting.  In  experiments  with  the  1.23  cm  tube,  the 

microphone  was  an  adapted  one-half  inch  B-K  condenser  microphone  with 

-3 

diaphragm  made  of  Mvlar  of  1.27  x  10  cm  thick.  Both  microphones  were 
powered  bv  a  B-K  tvpe  2801  power  suoplv  and  thpir  outputs  were  nreamnli- 
fied  bv  a  B-K  tvne  2815  cathode  follower  amplifier. 


The  speaker  used  in  experiments  on  stable  erases  was  one  previously 


used  bv  Bass  and  Detsch  in  their  investigation  on  the  sound  amnl i f icat i on 
in  gas  phase  reaction  svstems.^  The  polarizing  voltage  of  the  sneaker 
was  maintained  bv  a  Heathkit  model  IP-17  regulated  power  supnlv  and 
driven  bv  tone  bursts  generated  bv  a  GP  model  13PBB  tone  burst  generator 
and  amplified  by  a  Hewlett  Packard  (HP)  Model  467A  amplifier. 

Discharge  Circuit 

The  discharge  circuit  used  in  this  experiment  was  a  modification  of 

28 

the  BCR  discharge  circuit  designed  by  M.M.T.  Lov  and  P.A.  Roland.  The 
circuit  diagram  is  shown  in  Figure  2.  The  original  lf>  F  Mylar  capacitor 
was  replaced  by  three  ID  pF  Micamold  capacitors  set  in  parallel.  The 
nulse  transformer  was  an  automobile  ignition  coil,  Pelco  D523,  and  the 
SCR  a  tvoe  50  RIA100.  The  2  volt-lp  second  trigger  pulse  was  generated 
by  a  HP  model  214A  nulse  generator.  The  transformer  was  placed  near  the 
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sound  cube  and  connected  to  the  SCR  circuit,  housed  in  a  separate  box,  bv 
a  shielded  cable.  The  discharge  circuit  was  able  to  produce  a  discharge 
of  approximately  2.5  joules  in  1.23  cm  tube,  30  to  50  cm  long  at  pressure 
of  40  Torr. 

Electronics 

An  external lv-triggered  model  214a  HP  Pulse  generator  was  used  to 
trigger  the  discharge  and  the  Nicolet  oscilloscope  trace.  The  large 
capacitor, at  relatively  high  voltage,  then  discharged  through  the  primary 
of  the  transformer.  The  voltage  oulse  from  the  secondary  of  the 
transformer  produced  the  discharge  in  the  gas  sample  and  a  metastable 
state  with  an  overpopulation  of  excited  vibrational  states  resulted. 

Efforts  to  generate  the  sound  pulse  in  the  metastabale  gas  using  a 
speaker  in  the  end  of  the  discharge  tube  were  unsuccessful  due  to  the 
masking  effects  of  the  discharge  on  the  sound  wave  and  damage  done  to  the 
transducers  by  the  discharge.  These  difficulties  were  finally  solved  by 
placing  a  metal  screen  a  few  cm  from  the  speaker/reflector,  as  discussed 
in  the  preceding  section.  In  this  configuration,  the  sound  pulse  was 
generated  by  a  sudden  increase  in  temperature  of  the  gas  in  the  region 
where  the  discharge  occured.  The  increase  in  temperature  caused  a  pres¬ 
sure  increase  which  expanded  into  the  cool  section.  The  expansion 
produced  a  pulse  which  reflected  back  and  forth  in  the  tube.  The 
pressure  pulse  was  detected  each  time  it  was  reflected  from  the  micro¬ 
phone  end  of  the  tube.  The  signal  from  the  microphone,  after  being 
amplified,  was  passed  through  an  Ithaco  4302,  24  dB/octave  filter  and 
then  applied  to  the  vertical  input  of  the  Nicolet  digital  oscilloscope. 
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The  signal  thus  displayed  on  the  oscilloscope  screen  was  stored  in  the 
scope  memory  and  finallv  transferred  to  a  HP  model  85  computer  to  He 
analyzed . 

The  sampling  time  of  the  oscilloscope  may  be  adiusted  from  500  nsec 
to  200  sec  per  ooint.  In  this  experiment  the  sampling  time  was  set  at  10 
usee  per  point.  The  total  number  of  points  recorded  and  stored  bv  the 
oscilloscope  was  4096  points.  The  data  points  of  reflected  pulses  were 
transferred  to  the  computer  and  Fourier  analvzed  one  pulse  at  a  time. 

Gas  Handling  S vs  tern 

The  system  was  evacuated  with  a  Welch  Duo  seal  fore  pump  which  was 
connected  to  the  manifold  of  the  gas  handling  system.  The  pressure  of 
the  system  was  measured  with  M.K.S.  Baratron  type  170  M-25  and  PDR-D-1 
transducers.  The  first  being  used  for  the  pressure  range  up  to  10  Torr 
and  the  latter  used  from  10  to  100  Torr.  The  fore  nump  was  capable  of 
evacuating  the  system  to  less  than  2  millitorr.  Gas  samples  were  stored 
in  storage  tanks  and  connected  to  the  sound  tube  through  a  manifold  and  a 
series  of  valves. 

EXPERIMENTAL  PROCEDURE 

The  absorption  and  velocity  measurements  were  made  in  mire  V  and 
mixtures  of  with  He  and  H^.  The  gas  mixtures  were  made  using  the 
method  of  partial  pressure.  A  storage  tank  was  alternatelv  filled  with 
N9  and  He/H?  until  the  desired  concentration  at  total  pressure  of  one 
atmosphere  was  established.  The  gases  were  then  allowed  to  mix  at  least 
one  dav  before  measurement  was  made.  Since  the  relaxation  time  and  the 
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maximum  relaxation  frequency  of  N0  are  very  sensitive  to  the  concentra¬ 
tion  of  imnurities,  especially  water  vanor,  it  was  necessary  to  ensure 
the  purity  of  the  gas  under  consideration.  For  this  nuroose  the 
following  steps  were  followed.  Prior  to  making  any  measurement,  the 
svstem  was  thoroughly  evacuated  and  the  combined  outgassing  and  leaking 
rate  was  frequently  checked.  (This  svstem  was  able  to  be  evacuated  to 
about  2  milli  Torr  and  the  combined  outgassing  and  leaking  rate  was  less 
than  0.5  milli  Torr  per  minute  at  a  pressure  of  20  Torr.)  The  pressure 
in  all  pipes  connecting  the  storage  tanks  and  the  manifold  of  the  gas 
handling  svstem  was  always  kept  lower  than  the  pressure  of  interest  in 
the  sound  tube  to  avoid  impuritv  leaks  into  the  sound  tube.  The  sound 
tube  was  then  flushed  several  times  with  the  gas  sample.  In  doing  so, 
the  evacuated  sound  tube  was  filled  with  the  test  gas  to  about  100  Torr 
in  pressure  and  re-evacuated.  This  was  repeated  several  times.  The 
electronic  equipment,  except  the  discharge  circuit,  was  warmed  un  before 
measurements  were  made.  The  band-pass  filter  was  set  to  | ass  frequencies 
from  2  to  6.3  kHz. 

When  everything  was  readv,  the  pump  was  closed  off  from  the  svstem 
and  the  gas  sample  to  be  tested  was  introduced  into  the  tube  until  the 
desired  pressure  was  obtained.  After  thermal  equilibrium  with  the  tube 
wall  was  established,  the  pressure  and  the  tube  wall  temperature  were 
recorded  and  the  discharge  circuit  was  turned  on  with  the  capacitor 
voltage  set  at  a  desired  value.  To  avoid  impurities  venerated  bv  the 
sparks,  the  sound  signal  was  taken  from  the  first  spark  that  gave  tKe 


expected  signal  shane.  To  further  ensure  the  ouritv  of  the  gas  sarnie 
under  consideration,  every  sound  signal  recorded  was  for  a  fresh  gas 
sample  that  had  been  in  the  tube  no  longer  than  60  seconds.  tv,p 
procedure  discussed  above  was  repeated  for  measurements  at  different 
pressure  and  tube  lengths  (the  distance  between  reflectors').  ror  a  given 
set  of  these  parameters,  data  was  taken  for  approximately  nine  discharges 
and  averaged. 


FOHRIRR  TRANSFORM 


The  absorption  and  velocity  of  sound  at  a  given  freouencv  were 
determined  bv  Fourier  analyzing  the  sound  pulse  after  each  reflection. 
For  this  nuroose,  the  signal  detected  bv  the  microphone  was  filtered  and 
then  digitized  by  the  oscilloscope.  The  starting  time  and  ending  time 
for  each  of  the  reflected  pulses  were  read  visually  off  the  oscilloscope 
and  stored  in  the  computer.  The  digitized  signal  between  these  two  times 
was  then  transferred  to  the  computer  and  Fourier  analyzed.  The  fundamen¬ 
tal  freauencv  in  the  analysis  was  set  at  390.6  bv  setting  the  tine 
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interval  covered  in  each  analysis  at  256  points  or  2.56  x  in  seconds. 
If,  for  example,  there  were  N  digitized  points  between  the  beginning  and 
end  of  the  pulse  as  read  visually  from  the  oscilloscope,  the  signal 
Fourier  analyzed  would  consist  of  N  points  transferred  from  the  oscillo¬ 
scope  memory  followed  bv  256-N  zero  points.  The  sampling  period,  T  =  10 
psec,  easily  satisfied  the  condition  that  1/  T  =  2  x  maximum  freauencv 

since  the  highest  freauencv  with  an  appreciable  component  in  the  oulse 
29  30 

was  6  k"z.  ’  The  number  of  points,  256.  was  limited  bv  the  time 
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required  bv  Che  computer  to  make  the  analysis.  The  290  Hz  frequency 
spacing,  however,  was  quite  adequate  for  the  purpose  of  this  experiment. 

The  average  sound  velocity  for  one  round  trio  in  the  tube  was 
obtained  bv  the  equation 
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where  a  is  the  absorption  coefficient  (neper/cm)  and  6  is  the  effective 
reflection  coefficient.  Bv  making  measurments  at  more  than  one  value  of 
L,  both  a  and  6  can  be  determined.  Once  the  reflection  coefficient  is 
determined,  a  can  be  determined  from  the  amplitude  without  varying  the 
length  of  the  sound  oath. 

11 

A  computer  program  developed  by  the  Hewlett  Packard  Company  was 
modified  in  order  to  determine  the  sound  absorption  and  velocity  of  each 
frequency  component  of  interest  directlv  from  the  amplitude  and  phase 
obtained  from  the  Fourier  analvsis.  The  computer  program  was  developed 
bv  using  the  Fast  Fourier  Transform  (FPT)  algorithm  that  tranformed  2'’ 
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sample  points  by  means  of  an  N  point  transform.  Tbe  results  of  this 
procedure  were  checked  bv  direct  measurements  of  the  sound  absorption  and 
velocitv  of  the  pulse  after  it  passed  through  the  filter  set  at  a  single 
frequency.  The  two  methods  gave  results  that  agreed  within  experimental 
error.  The  computer  program  used  to  Fourier  analvg.e  tb-e  sound  signal  is 
presented  in  the  appendix. 
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RESULTS  OF  MEASUREMENTS  IN  STABLE  GASES 

Sound  amplification  and  translational  temperature  were  obtained  bv 
comparing  the  measured  absorotion  and  velocity  in  excited  ?ases  with  the 
absorption  and  velocity  calculated  for  the  unexcited  gases  using  the 
numerical  solution  discussed  above.  To  check  the  accuracy  of  t^ese 
calculations  the  sound  absorption  and  velocity  were  measured  in  unexcited 
nitrogen  gas  and  mixtures  of  25%  Me  -  75%  N?,  50%  He  -  50%  N?  of  various 
pressures  and  frequencies.  In  these  measurements,  the  gas  was  contained 
in  a  cylindrical  glass  tube  7. 54  cm  inside  diameter.  A  tone  hurst, 
produced  hv  the  speaker  in  the  end  of  the  tube  was  allowed  to  propagate 
back  and  forth  between  the  speaker  and  the  reflector  while  the  amplitude 
was  monitored  bv  the  microphone.  The  amplitude  of  the  tone  burst  after 
it  was  reflected  m  times  is  given  bv 


A  =  A  e 
o 


-(2La+3  )m 


(22) 


where  a 

,  3  and  L  have 

been  defined 

previously. 
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corresponding  Co  a,  and  they  are  diso1 seed  from  each  other  bv  an  amount 
corresponding  to  s.  The  results  obtained  from  a  least  square  fit  of  a 
straight  line  to  the  data  are  given  in  Tables  3  and  4  along  with  the 
calculated  values.  From  comparison  of  values  in  these  tables  it  is 
concluded  that  individual  measured  and  calculated  values  of  absorption 
differ  bv  at  most  -2’-.  Note  that  relaxation  absorntion  in  nitrogen  »as 


at  room 

temperature  is 

ve  rv  sma 1 1 . 
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not 
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The  velocitv  was  obtained  bv  measuring  the  time  the  sound  took  to 
travel  a  distance  of  2L  (one  round  trip)  in  the  tube.  The  results  are 
compared  with  theoretical  in  Tables  5  and  A.  From  comparison  of  these 
values  it  is  concluded  that  the  maximum  difference  in  individual  measured 
and  calculated  values  of  velocities  is  -0.3". 

SOUND  ABSORPTION  AND  VF.UOCTTY  IN  FxdTF.n  MF.T  A  ST  ARLF  STATF.  OASFS 

Absorption  and  velocitv  measurements  in  excited  netastahle  state 
gases  were  made  with  two  different  tube  sixes,  2.34  cm  and  1.23  cm  inside 
diameter.  Cases  used  with  the  bigger  tube  were  pure  N0  and  mixtures  of 
2 ST  He  -  7  3"  N0  and  30"  He  -  30"  all  at  pressures  of  20  Torr.  In 
experiments  with  the  smaller  tube,  the  above  erases  and  a  10T  1^  -  POT  N,, 
fixture  at  pressures  of  20  and  40  Tirr  were  used.  Fxnerinents  were  begun 
hv  discharging  the  canacitor  through  the  orimarv  of  the  induction  coil  , 
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TABLE  5 


Sound  velocity 

Pressure  = 

in  stable  at  various  frequencies 

20  Torr,  tube  radius  =  1.27  cm. 

f (kHz ) 

V  (m/s) 

obs. 

V  (m/s) 

calc. 

1.6 

— 

342.4 

I 

2.0 

342. 8 

343.  1 

3. 15 

345.2 

344.3 

4.0 

344.8 

344. « 

B 

TABLE  6 

■ 

Sound  velocitv  in  stable 

50%  He  -  50%  N2 

mixture  at  various 

f 

=  2  kHz,  tube  radius  = 

1.27  cm. 

P(Torr) 

V  (m/s) 

obs. 

Vcalc.(m/S) 

20 

468.5 

469.7 

30 

471.6 

471.8 

40 

472.4 

473.  1 
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thus  producing  the  discharge  in  the  gas  sample.  The  electrical  discharge 
rapidly  increased  the  gas  pressure  in  the  region  between  the  two  elec¬ 
trodes  (the  region  between  the  reflector  and  the  screen;  see  Figure  1). 
The  rapid  increase  in  the  pressure  in  the  discharge  region  creates  a 
pressure  pulse  at  the  hot-cold  iunction  where  the  discharge  terminated  in 
the  gas.  This  pulse  travelled  back  and  forth  in  the  tube  and  was  moni¬ 
tored  bv  the  microphone.  The  microphone  response  to  the  sudden  pressure 
increase  is  illustrated  in  Figure  4.  In  this  figure,  positive  pressure 
corresponds  to  a  negative  voltage.  The  spikes  represent  the  sound  pulse 
that  reflected  back  and  forth  in  the  tube.  The  shape  of  the  response 
curve  without  the  superimposed  pulses  is  the  response  of  the  microphone 
to  the  rapid  increase  in  Pressure  that  decavs  with  a  decav  time  eaual  to 
that  for  cooling  the  gas  in  the  tube,  bv  filtering  and  Fourier  analyzing 
the  signal,  the  amplitude  and  the  phase  of  the  different  freouencv  compo¬ 
nents  in  the  pulse  as  a  function  of  time  can  be  determined.  Examples  of 
a  filtered  signal  are  shown  in  Figures  5  and  6.  The  legend  for  these  and 
other  figures  specifies  the  energy  in  the  discharge  in  joule/mole.  The 
energy  was  obtained  by  dividing  the  energv  stored  in  the  capacitor  before 
discharge  by  the  moles  of  the  gas  heated  bv  the  discharge.  It  is  esti¬ 
mated  that  approximate lv  70%  of  this  energv  actuallv  went  into  bating 
the  vibrational  modes  of  the  gas. 

The  decreasing  peaks  in  Figure  5  represent  the  sound  nulse  detected 
after  it  travelled  successive  round  trips  down  the  tube  (a  distance  of 
2  1).  The  first  peak,  is  the  original  pulse  detected  after  it  travelled  a 
distance  of  aoproximatelv  4  cm  through  the  cold  region.  The  second  peak 
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is  the  pu!se  detected  after  it  travelled  one  round  trip.  The  third  peak 
is  the  pulse  detected  after  it  travelled  two  round  trios,  etc.  The 
distance  between  two  adiacent  peaks  corresponds  to  the  tine  the  nulse 
took  to  travel  a  distance  of  2  L.  In  general  the  first  pulse  considered 
in  our  analvsis  was  the  one  reaching  the  microphone  after  travelline  one 
round  trip  in  the  tube. 

The  average  absorption  per  unit  leneth  between  any  two  successive 
reflections  is  given  bv 


a  =  ( -In — r - 3)/2L;  m  =  1,2,3,..., 


The  absorption  determined  bv  Eauation  23  was  assumed  to  be  the  absorption 
at  a  time  equal  to 


t  =  (t  +t  ) / 2 , 
m  m+1 


where  t  is  the  time  after  the  nulse  was  reflected  m,  times,  taken  at  the 
m 

center  of  the  pulse.  In  the  case  that  the  sound  pulse  is  amplified,  a  is 
time  dependent  because  the  vibrational  temperature  is  time  dependent.  In 
experiments  with  the  2.54  cm  diamter  tube,  B  was  obtained  from  the  exper¬ 
iment  on  stable  gases  as  discussed  in  the  preceding  section.  Mith  a  1.23 
cm  diameter  tube,  attempts  to  determine  3  bv  using  a  tone  burst  from  the 
speaker  failed  because  the  sound  attenuation  at  low  pressure  was  too 
large.  The  tone  burst  decayed  out  verv  ranidlv.  In  this  case,  a  sound 
pulse  generated  by  the  discharge  was  used  instead  of  the  tone  burst.  3 
was  determined  from  the  intercept  of  the  plot  of  the  logarithm  of  amplitude 
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versus  distance,  L,  as  discussed  previously.  In  this  case,  of  course, 
the  discharge  circuit  was  adjusted  so  that  the  energy  ner  unit  volume  in 
the  discharge  remained  constant  for  all  reflector  senarations. 
Discharges  were  chosen  for  different  reflector  distances,  L,  so  that  the 
amplitudes  for  the  first  reflected  pulse  were  the  same.  Table  7  shows 
the  coefficient  of  reflection  in  Neper /ref lection  averaged  over  fre¬ 
quency  components  from  1.9  to  4.2  kHz,  observed  in  gases  used  in  this 
experiment . 

TABLE  7 

Coefficient  of  reflection  (8)  observed  in 
experiments  with  1.23  cm  diameter  tube. 


(Neper/ ref lection) 


GAS 

P  =  20  Torr 

P  =  40  Torr 

N2 

0. 1 6±0. 07 

0.15+0.03 

25He-73N9 

0.05+0.03 

0.09±0.02 

50He-50N2 

0.10+0.04 

O 

• 

1+ 

o 

• 

o 

10H2-90N2 

0.03+0.01 

c 

9 

o 

+1 

r- 

o 

c 

Results  of  Measurements  in  Excited  Metastable  Pure  Nitrogen 

Sound  absorption  and  velocity  observed  in  excited  metastable  nure 
nitrogen  are  nreeented  in  Tables  8-13.  Table  8  gives  the  average 
absorption  coefficient  over  the  total  time  the  pulse  was  observed.  Also 


TABLE  1  1 


mencs  it  is  estimated  that  the  measured  values  stiver!  in  Table  8  are  vood 
to  -3%.  Aa  is  the  difference  between  the  observed  and  calculated  absorp¬ 


tion  ( Aa  =  a  ,  -a  .  ). 

calc.  obs. 

It* 

T^e  eain  ( Aa)  in  dB/cm  should  be  frequencv  indenendent.  Making;  use 
of  this  fact,  bv  averaging  Aa  over  all  frequencv  components  in  Tabl e  B, 
we  found  Aa  of  approximately  0.003  dB/cm  in  the  bigger  tube,  and  0.nQ4 
dR/cm  and  0.005  dB/cm  in  the  smaller  tube  at  pressures  of  20  and  40  Torr 
respectively.  This  difference  is  only  slightlv  more  than  experimental 
error.  If  it  is  due  to  relaxation,  it  should  be  changing  with  time. 
Therefore,  in  Tables  9-13  the  absorption  and  velocity  between  individual 
reflections  is  tabulated  as  a  function  of  the  time  at  the  midpoint  of  the 
interval  considered.  Any  change  in  velocity  and  absorption  with  time  was 
within  experimental  error. 

At  first  sight  it  might  be  somewhat  surprising  that  the  temperature 
of  the  nitrogen  gas  was  not  changed  more  bv  the  discharge.  (The  velocity 
measurements  show  the  gas  temperature  to  be  only  a  few  degrees  above  the 
tube  wall  (room)  temperature.)  In  later  experiments  with  mixtures,  the 
discharge  was  found  to  heat  the  gas  from  10  to  40°C.  However,  a  closer 
look  at  the  properties  of  pure  nitrogen  explain  the  absence  of  the 


heating.  'Mahan's  calculations  from  measured  electron  cross  sections 


in  nitrogen  indicate  that  for  the  voltages  and  nressures  of  our 
discharge,  practically  all  of  the  discharge  energy  goes  into  vibration. 
Further,  the  vibrational  relaxation  time  in  nure  nitrogen  is  so  much 
greater  than  the  thermal  diffusion  time  in  the  gas  that  the  vibrational 
energy  is  conducted  to  the  walls  without  appreciably  heating  the  gas. 

Though  the  difference  between  measured  and  calculated  absorption  is 

small,  the  observed  values  were  consistantlv  smaller  than  calculated  as 

would  be  expected  if  amplification  were  present.  If  this  difference 

between  measured  and  calculated  absorption  is  due  to  relaxation  it  is 

10 

larger  than  Shields'  calculation.  Assuming  the  vibrational  temperature 
is  between  1000  and  2000°K  during  the  measurement,  the  gain  would  indi¬ 
cate  a  relaxation  time  in  nitrogen  of  the  order  of  a  few  milliseconds  per 
atmosphere.  This  relaxation  time  is  considerably  shorter  than  the 
relaxation  time  in  pure  nitrogen  and,  for  the  estimated  amount  of  energy 
in  the  discharge,  the  vibrational  temperature  could  not  exceed  S0no°K. 

The  small  discrepancy  is  thus  unexplained,  and  is  assumed  to  be  due 
to  a  consistent  experimental  error.  In  experiments  on  mixtures  of  nitro¬ 
gen  with  helium  and  hydrogen,  which  will  be  discussed  later,  Act  here 
observed  will  be  used  as  a  correction  factor. 

The  uncertainty  in  velocity  at  a  particular  time,  estimated  from  the 
standard  deviation  in  the  individual  measurements,  is  approximately 
-0.4TS.  Tn  what  follows,  the  temperature  of  the  gas  will  be  obtained  hv 
comparing  measured  velocities  with  values  calcnlated  as  functions  of 


temnerature  by  numerical  solutions.  Tomneratures  obtained  in  this  wav 
are  estimated  to  be  accurate  to  -1*C. 

Results  of  Measurements  i n  He-N'.,  Mixtures 

Sound  absorption  measurements  in  He-N^  mixtures  are  presented  in 
Tables  14-17.  The  variation  in  Aa  with  time  in  50%  He  -  50%  N  mixtures 
is  plotted  in  Figure  7. 

The  addition  of  He  molecules  to  gas  drastically  reduces  the  relaxa¬ 
tion  time  of  N? .  In  a  25%  He  -  75%  N  mixture  th»  relaxation  time  is 
reduced  to  about  1040  and  520  milliseconds  at  pressures  of  20  and  4!)  Torr 
respectively,  and  in  a  50%  He  -  50%  \  mixture  to  520  and  260 
milliseconds.1*’  However,  compared  to  the  thermal  conduction  relaxation 

time  (t  )  which  is  of  the  order  of  in  milliseconds,  the  ratio  of  x  and 
tc  v 

x  is  still  large.  The  vibrational  energy  transferred  into  translation 
tc 

during  the  transit  time  is  still  small,  and  the  gain  Aa  is  expected  to  be 
so  small  that  it  will  be  difficult  to  observe.  The  observed  value  of  Aa 
given  in  Tables  14  and  15  was  corrected  by  subtracting  the  Aa  observed  in 
pure  nitrogen.  When  this  is  done,  no  appreciable  gain  is  observed  in  the 
25%  He  -  75%  mixtures,  nor  is  there  an  observable  change  in  it  with 
time. 

However,  in  the  50%  He  -  50%  N?  mixtures  there  is  a  measurable  gain. 
Averaged  overall  frequency  components  in  Table  16,  Aa,  after  correcting 
bv  subtracting  Aa  observed  in  N  ,  is  0.0018  dB/cm  and  0.0074  dB/cm  at 
pressures  of  20  and  40  Torr  respectively.  The  measured  Aa  values  are 


Time  average  sound  absorption  (a)  in  excited  metastable  25%  He  -  75%  N„  mixture 
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7.  Gain  (Act)  as  a  function  of  time  observed  in  50%  He  -  50%  V 
at  pressure  40  Torr;  tube  radius  =  0.615  cm;  discharge  *ner<™ 

ioule/mole.  Rolid  line  is  calculated  gain  setting 
)/dT  =  0.05,  T  =  (3000 0 V  at  4  msec;  x  =  7.7  msec. 


m  -j 


53 


plotted  versus  time  in  Figure  7.  The  magnitude  and  time  dependence  of 
this  gain  is  comnared  to  theorv  below. 

Tables  13  and  Id  give  renresentative  velocities  measured  in  25%  ”e  - 
75%  and  50%  He  -  50%  mixtures  at  a  pressure  of  40  Torr  and  tube 

length  of  30  and  40  cm  resnect  i  velv.  This  time  the  measured  velocities 
show  an  appreciable  increase  in  translational  temperature.  This  indi¬ 
cates  that,  with  the  added  He,  a  larger  amount  of  the  discharge  energy  is 
going  directly  into  the  translational  energy  of  the  gas  mixture.  The 

logarithm  of  the  difference  between  the  translat ional  temperature  (T  )  as 

r 

determined  from  the  velocity  measurements  and  the  tube  wall  temperature 

(T  )  is  plotted  versus  time  in  Figures  3  and  d.  The  decav  times  calcu- 

lated  from  the  curves  are  compared  with  theoretical  values  of  x  in 

tc 

Table  20.  The  measured  values  of  x  do  have  roughlv  the  expected 

variation  with  pressure  and  tube  radius  as  predicted  bv  Eauation  Id. 

In  mixtures  in  which  x  was  considerably  longer  than  the  time  for 

tc 

the  sound  to  make  one  round  trio  in  the  tube  (25%  He  -  75%  and  in? 

-  d0%  at  a  pressure  of  40  Torr)  the  first  data  point  of  the  transla¬ 
tional  temperature  (see  Figures  8  and  12)  vras  far  off  the  straight  line. 
This  was  attributed  to  a  non-uniform  heating  of  the  gas  by  the  discharge. 

When  x  was  long,  compared  to  the  transit  time  for  the  sound,  there  was 
tc 

not  time  for  thermal  diffusion  to  smooth  out  the  hot  and  cold  regions  in 
the  gas  left  bv  the  non-uniform  heating  of  the  discharge.  The  hot  snots 
evidently  raised  the  measured  sound  velocity.  At  lower  pressure  the 
effect  was  not  so  obvious  due  both  to  the  fact  that  the  discharge  was 
more  uniform  and  the  thermal  diffusion  time  was  shorter. 


TABLE  19 


Sound  velocity  in  m/s  as  a  function  of  time  in  excited  metastable 

50%  He  -  50%  N?  mixture  at  Pressure  =  40  Torr;  tube  length  =  40  cm; 

3 

tube  radius  =  0.615  cm;  discharge  energy  =  16.7  x  10  ioule/mole. 


Time 
\ms  ) 
f  \ 

(Hz)  \ 

1 

1 

2.91 

{ 

4.53 

i 

- 1 

6.21 

1 

1  V  at 

C  <3  1  C. 

T=303°K 

1953 

1  1 
471.8  | 

I 

468.4 

467.2 

| 

467.6 

2344 

t 

471.5  | 

468.6 

r  1 

467.5 

469.2 

2734 

471.3  | 

I 

468.9 

468.3 

470.4 

3125 

[  ! 
471.3  | 

| 

469.9 

469.3 

471.3 

3515 

1 

471.7  | 

- h 

471.1  | 

1 - ^ 

472.2 

3906 


472.  1 


471.2 


470.4 


472.° 


Time  (msec) 


P"[Cin?F,  8.  The  difference  between  translational  and  wall  temperature  (T 
-  T  )  as  a  function  of  time  observed  in  25°  ue  -  7  5"  mixture.  ^ 

pressure  (°)  =  20  Torr;  tube  rad’us  (21  =  0.015  cm;  tube  lenfTtv'  0.1  =  i'1 
cm.  v  -  n  =  20  Torr;  p  =  0 . 5 1  5  c  n ;  I .  =  4  3  c  m .  2  —  e=  40  Torr;  P  = 

n .  b  i  i  cm;  i.  =  30  c"1 .  a  —  p  =  in  Torr;  2  =  0 .  is  l  5  cm ;  1  —  03  cm.  c  —  n  = 

20  Torr;  D  =  1.27  cm;  1.  =  50  on.  Solid  lines  are  least  senates  fit  to 
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FTW’F  9.  The  difference  between  translational  and  wal 
-  T  )  ns  a  function  of  time  observed  in  Sf'7  He  -  5*'1"  V 
2n  f o r r ;  P  =  d.AlS  cm:  !.  =  3°  cm.  □  -  D  =  in  Torr; 
cm.  o  —  P  =  2<3  torr:  p  =  1 .  ?7  cm;  i.  =  am  cm.  c.ii  i 
snuares  fit  to  data. 
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TABLE  20 

Measured  values  of  the  decay  time  for  translational  temperature  in 
0.615  cm  radius  tube  compared  with  values  of  x  calculated 
by  Equation  19.  Values  are  given  in  milliseconds. 
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-ry 
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Obse  rveH 

Calculated 

■ 
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5.78 
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25%  He  -  75%  N2 
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Results  of  Measurement  In  H^-Mn  Mixtures 

In  10%  -  00’/  N^  mixtures  the  vibrational  relaxation  times  are 

reduced  to  about  106  and  53  milliseconds  at  pressures  of  20  and  40  Torr 
16 

respectively.  The  thermal  conduction  relaxation  times  in  this  mixture 

at  these  pressures  are  approximately  5  and  10  milliseconds  for  the  1.23 

cm  tube  diameter  as  predicted  hv  Foliation  19.  Thus  the  ratio  x  /x  is 

v  tc 

relatively  small.  For  a  vibrational  temperature  of  3000°K,  the  gain,  Aa, 
is  predicted  by  Shields'  calculation  to  be  about  0.009  HR.Vn  at  a  pres¬ 
sure  of  40  Torr. 

The  observed  time  averaged  absorption  and  main  (after  subtracting  bv 
Aa  observed  in  pure  nitrogen)  are  presented  in  Table  21.  As  expected 
from  Shields'  calculations  the  gain  is  independent  of  frenuencv  to  the 
accuracy  of  measurements.  Table  22  and  23  give  the  observed  velocity  at 
20  and  40  Torr  pressures.  The  corresponding  translational  temperatures 
are  plotted  in  Figure  11.  As  shown  in  Table  21,  the  average  gain  was 

found  to  be  0.02  dB/cm  at  20  Torr  and  0.017  dB/cm  at  40  Torr.  Figure  1° 

shows  the  variations  in  Aa  versus  time  at  40  Torr.  At  20  Torr,  the  time 
dependent  value  of  Aa  showed  too  much  scatter  to  allow  a  reliable  deter¬ 
mination  of  the  decav  time.  This  was  attributed  to  the  fact  that  x  was 

vc 

shorter  than  the  sound  transit  time  over  several  reflections. 

COMPARISON  OF  THFORFTICAL  AND  MFASl!RFP  PPLAX ATI OV  Tl'TS  fv  tiif  1  or  _ 

90%  n2  mixturf. 

When  x  and  x  are  of  comparable  signs  tv,e  iecav  tine  of  the 
v  tc 

translational  temperature  should  he  longer  than  t^er-il  diffusion  tine 


60 


TABLE  21 

Tinie  average  sound  absorption  (a)  and  gain  (Aa  =  a  ,  -  a  ,  )  in  excited 

calc  obs 

metastable  10%  H  -  90%  N?  mixture.  Tube  radius  =  0.615  cm; 

3 

discharge  energy  =  44.8  x  10  ioule/mole  at  20  Torr  and  22.4  x 

3 

10  ioule/mole  at  40  Torr.  Aa  has  been  corrected  bv 
observed  in  pure  N_.  Values  are  given  in  dB/cm. 


TABLE  22 


Sound  velocity  in  m/s  as  a  function  of  time  in  excited 
10%  -  90%  at  Pressure  =  20  Torr;  tube  length 

3 

tube  radius  =  0.615  cm;  discharge  energv  =  48.8  x  10 
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TABLE  23 


Time  (msec) 


FIGURE  10.  Gain  (ia)  as  a  function  of  time  observed  in  1 0J<  -  Rd”/ 

mixture  a^  pressure  =  40  Torr;  tube  radius  =  0.615  cm;  discharge  ener^v  = 
48. 8  x  10  ioule/mole.  Oata  noints,  O-  observed  at  tube  length  =  30  cm; 
A  -  observed  at  tube  length  =  43  cm.  Curved  is  calculated  gain  setting 
d(lnr  )/dT  =  0.05.  Curve  3  is  4T  x  Id  computed  from  Filiation  37. 
CurveVC  is  the  calculated  gain  bv  setting  d(lnx  ~^)/dT  =  O.n;. 
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FIGURE  11.  The  difference  between  trnnslation.nl  and  wall  temperature 
(AT)  as  a  function  of  time  for  experiment  on  10K  H  -  80?<  \'0  mixture  at 
pressure  =  20  Torr;  tube  radius  =  n.315  cm;  discharge  eneriv“=  48.3  •<  103 
joule/mole.  Oata  points,  O  -  observed  at  tube  length  =  40  cm;  A  - 
observed  at  tube  length  =  43  cm.  Solid  curve  is  the  theoretical  curve 
computed  from  Equation  34. 


due  to  the  feeding  of  energy  into  translation  from  vibration 


The 


difference  is  noted  between  the  measured  and  calculated  relaxation  times 
for  10%  -  00%  N  mixtures  in  Table  20.  The  energy  truing  from  vibra¬ 
tion  to  translation  also  makes  the  decay  time  for  vibrational  temnerature 

shorter  than  the  calculated  t  .  The  relationship  between  x  as  calcu- 

tc  tc 

lated  from  Equation  19  and  these  decay  times  for  translational  and 
vibrational  temnerature  as  determined  from  measured  velocities  can  be 


obtained  as  follows.  Assume  the  translational  mode  gains  energy  from  the 
vibrational  mode  and  loses  it  to  the  wall.  The  time  rate  of  change  of 
the  translational  and  vibrational  energy  can  be  written  as 
dE  C . (T  )  C 

— tr  =  — v_  ( T  _T) - (T-T  ),  (25) 

dt  xv  v  xtc 


dE  C . (T  )  C . (T  ) 

— -  = - - — —  (  t  -T'l - i — —  (T  -x  ),  (26) 

dt  tv  i  v  w 

v  vc 

where  E  is  the  translational  energy,  E  is  the  vibrational  energv,  t 
tr  v  v 

T,  and  T  are  the  vibrational,  translational,  and  wall  temnerature 
w 

respectively,  is  the  vibrational  snecific  heat,  and  is  the  time 

independent  specific  heat  enual  to  6/2  R  per  mole  for  N.  and  H„.  x  ,  x 

2  2  v  vi¬ 


and  x  are  the  same  as  defined  before.  ev  setting 
tc 


T  -T  =  (T  -T  )  -  (T-T  ), 


dr. 
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tube 


Time  (msec) 


FIOURF  12.  The  difference  hetweer.  translational  and  wall  temnerature 
(AT)  as  a  function  of  time  for  experiment  on  10”  o  ,  -  bO7  \T o  mixture  at 
oressure  =  40  Torr;  tube  radius  =  0.61  ^  cm;  discharge  eners'v  =  24.4  x  I'M 
ioule/mole.  Oata  ooints,  a  -  observed  at  tuKe  li'i-'t'-'  =  25  o"’;  o  - 
observed  at  tube  length  =  40  cm.  Folid  curve  is  t'ne  theoretical  curve 
computed  from  Equation  19. 


P 
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A„  and  A,  can  be  related  to  A,  and  A.  bv  substituting  Eauations  33 
3  4  12 

and  34  into  Eauations  25  and  2h.  Separating  the  tent's  containing 
exp(-Ajt)  and  those  containing  exn(-A^t)  one  obtains 


(35) 


and 


(3b) 


The  procedure  in  analyzing  the  10^  -  90%  data  then  is  as 

fol lows 

1.  Assume  an  average  value  for  the  vibrational  temperature,  T  ,  over 

24 

the  period  of  measurement,  and  determine  C.  from  C  using  the  eauation 

l  v 

C.  =  C  -C  . 

1  V  oo 

2.  Using  known  values  of  t  ,  and  t  (=t  ),  and  the  assumed  C. 

v’  tc  vc  i 

determine  A^  and  A^  from  Eauations  31  and  32. 

3.  From  the  values  and  slopes  of  curves  in  Figures  11  and  12  at 

f  =  1  ,  determine  A,  and  A„ ,  4  and  A.  as  explained  above, 

t  c  1  2  3  4 

4.  Calculate  AT  using  Eauation  34  over  the  measured  time  interval 

v 

and  compare  with  the  average  value  assumed  in  Step  1. 

5.  Repeat  until  the  values  of  AT  obtained  in  Step  4  are  consistent 

v 

with  values  of  T  assumed  in  Step  1. 
v 


Using  a 
AT  =  24’, 


value 

and 


of  C .  =  0.21R  corresponding  to  T  = 
i  v 

f-^-(  1  n/.T)  1  _1  = 


1S00°’<;  values  of 


9.05  milliseconds  correspond i ng  ro  value® 
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from  Figure  11  at  C  =  6  milliseconds;  calculating  xt  from  Fquation  19; 

setting  x  =  x  ;  and  obtaining  x  from  Reference  16  we  obtained 
vc  tc  v 


AT(°C) 


532. 5e 


-o. 1 79 t 


492e 


-n. i90t 


(37) 


and 


AT  CC)  =  532.5e_n*  179t  +  4731. 4e  n*190t 


( 3S> 


v 

where  t  is  in  milliseconds.  The  curve  in  Figure  11  drawn  through  the 
noints  taken  at  20  Torr  is  calculated  from  F,a nation  37.  In  a  similar  wav 
the  equations 


AT( °C) 


238. Oe 


-n.089t 


220. 3e 


-0. 1 lOt 


(39) 


and 


AT  (*C) 
v 


238. Oe 


— 0.089 1 


+ 


2622. 4e 


-0. 110t 


(40) 


fit  data  points  ohserved  at  40  Torr  pressure  in  Figure  12.  The  constants 

in  the  40  Torr  equation  were  obtained  bv  using  values  of  AT  =  24°C  and 

d(lnAT)  1  _  msec  as  determined  from  Figure  12  at  t  =  10  msec  and  bv 

dt 

setting  Ck=0.21R  corresponding  to  Tv=1500°K.  It  should  he  noted  that 
these  equations  have  been  obtained  from  the  velocitv  measurement  without 
reference  to  absorption. 

The  agreement  between  the  measured  and  calculated  time  dependence  of 
AT  in  Figure  11  is  impressive  since  the  onlv  adiustable  parameter  is  the 
vibrational  temperature. 
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In  Figure  12  the  time  dependence  is  correctly  predicted  after  approx¬ 
imately  10  msec,  but  the  velocitv  obtained  from  the  first  reflected  pulse 
is  too  high.  This  was  also  the  case  for  the  50%  He  -  50%  \'9  mixture  at 
40  Torr.  Here,  as  there,  the  difference  can  be  attributed  to  the  fact 

that  the  first  reflection  is  occurring  for  times  shorter  than  x  and  the 

tc 

temperature  variations  left  bv  the  discharge  have  not  been  smoothed  out 

nor  the  metastable  state  developed. 

F.quations  38  and  40  give  the  vibrational  temperature,  AT  ,  as  a 

function  of  time.  Hote  these  equations  have  been  determined  iust  from 

measured  AT  vs  time  values  and  known  values  of  x  ,  x  and  i  .  T^e 

v  tc  vc 

value  of  AT  predicted  at  t  =  x  is  2100I>K  for  20  Torr  and  1200°K  for  40 
v  tc 

Torr.  These  values  would  imply  that  S.4%  of  the  discharge  energv  was  in 

the  gas  at  that  time.  This  percentage  is  unreasonably  small.  The  gain 

measurements  in  the  next  section  indicate  a  higher  vibrational 

temperature  was  present.  This  discrepancy  would  infer  that  the  value  of 

t  in  H  -N„  mixtures  used  in  calculating  the  constants  in  Equation  _,q  and 
v  2  2 

40  which  were  obtained  from  fluorescent  measurements  by  Frev,  et  al.' 
were  too  small. 

COMPARISON  OF  THEORETICAL  AND  MEASURED  GAIN  IN  10%  H0  -  DOT  N  AND  SOT 
He  -  50%  N2.  " 

In  two  cases  it  has  been  possible  to  plot  the  vain  versus  time  (50% 

He  -  50%  N2  at  40  torr  and  10%  -  90%  N  at  40  Torr).  For  the  H0-\T9 

mixture  it  has  been  possible  to  get  T  vs  time  from  the  velocitv  mea'ure- 

v 

ments  as  well.  It  is  now  interesting  to  compare  these  measured  gains 
with  values  calculated  from  Shields'  theory.  Por  10%  H9  -  Q0o/  (  the 


measured  values  of  gain  vs  time  are  plotted  in  Figure  10  also  shown 

plotted  in  this  figure  to  a  different  scale  is  IT  vs  time  calculated  bv 

v 

Equation  40  (curve  3).  Since  both  T  and  T  are  determined  bv  Equations 

v 

39  and  40  as  a  function  of  time,  it  is  also  possible  to  calculate 

expected  gain  versus  time  using  Shields'  theorv.  The  result  is  plotted 

in  Figure  10  as  curve  C  assuming  value  -pr(lnx  ^ )  =  0.01.  The  result 

dT  v 

gives  magnitudes  of  gain  that  are  too  small  bv  a  factor  of  10  and  a 

negative  slope  that  is  too  large  bv  a  factor  of  two.  There  are  three 

parameters  that  affect  the  calculated  gain:  i  ,  -^-(lnx  *),  and  T  .  The 

v  dT  v  v 

magnitude  of  the  gain  can  be  increased  bv  increasing  -pr(lnx  ')  or  T  , or 

dT  v  v 

decreasing  x  •  However  the  slope  of  the  gain  with  time  is  changed  little 
v 

by  changing  anv  of  these.  The  measured  small  decrease  in  gain  with  time 

implies  a  T  that  changes  more  slowlv  with  time  than  predicted  bv  Eaua- 
v 

tion  40.  This  descrepancv  may  be  due  to  the  fact  that  most  of  the 

measurements  in  this  mixture  were  made  for  times  less  than  x  while 

tc 

the  temperature  T^  and  T  were  still  probably  far  from  uniform  in  the 
tube.  The  appropriate  magnitude  of  the  gain  would  be  predicted  bv 
Shields'  calculation  if  70%  of  the  discharge  were  in  vibration  at  t  =  x 

tc 

msec  and  -r=( lnx  ^ )  =  0.05.  Curve  A  in  Figure  1C  corresponds  to  values 
dT  v 

calculated  bv  Shields'  theorv  assuming  -pr(  lnx  =  0.05  and  T  equal  to 

dT  v  v 

3000°K  at  t  =  7.5  msec  which  decays  with  decav  tine  eaual  to  x  .  This 

tc 

would  correspond  to  70%  of  the  discharge  energy  residing  in  the 
vibrational  energy  of  the  gas  at  t  =  x  . 
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For  Che  case  of  50%  He  -  50%  mixture  ac  40  Torr,  it  was  not 

possible  to  get  T  vs  time  from  the  velocity  measurements.  However, 

v 

because  t  is  so  much  larger  than  t  in  this  gas  it  is  safe  to  assume 
v  tc 

T  decays  with  a  relaxation  time  equal  to  i  .  Assuming  the  sane  oer- 
v  tc 

centage  of  discharge  energy  goes  into  vibration  in  this  mixture  as 

with  the  H„  mixture,  a  value  of  T  at  t  =  i  would  eaual  6h00°K.  Using 
2  v  tc 

this  value  of  T  and  again  assuming  -4glnvT  *)  =  0.05.  The  calculated 
v  dT  v 

gain  in  50%  He  =  50%  N0  is  plotted  as  a  solid  line  in  Figure  7. 

AN  EXAMINATION’  OF  THE  CONSISTENCY  OF  THE  RESULTS. 

It  is  interesting  to  note  a  number  of  consistencies  in  the  data. 
First  consider  the  AT  vs  time  data  that  has  been  obtained  from  the 
velocity  measurements. 

1.  The  initial  translational  temperatures  that  result  from  the 

discharge  are  reasonable.  In  pure  N9  the  translational  temperature 

changed  onlv  slightly  as  expected  due  to  the  small  amount  of  discharge 

energy  going  directly  into  translation  and  the  big  value  of  t  /t  .  In 

v  tc 

the  mixtures,  the  value  of  AT^  increase  approximately  linearly  with  the 
discharge  energy  as  might  be  expected  (see  Figure  13"). 

2.  The  time  dependence  of  AT  in  all  mixtures  agrees  roughlv  with 

values  calculated  from  values  of  t  ,  r  and  i  obtained  independently 

tc  vc  v 

(Table  20). 

3.  The  gain  per  unit  length  was  found  to  be  independent  of  freauencv 
as  predicted  theoretically  for  the  freauencv  range  considered  (see  Table 


21  ) 


Discharge  Energy  x  10  (joule/mole) 


FIGURE  l?.  The  increase  in  translational  temperature  us  u  function 
of  discharge  energv  in  gases  under  studv.  Solid  line  is  visual  straight 
line  fit  to  data. 
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4.  In  both  10%  -  40%  and  50%  He  -  50%  ,  the  train  vs  time 

d  ,  -1 

curve  indicates  the  tenoerature  deoendence  of  the  transition  rate  * 

is  given  bv  -^rln(x  S  =  0.05  and  70%  of  discharge  energv  is  in  the 
dT  v 

vibrational  mode  of  the  gas. 

There  are  two  aspects  of  the  data  still  unexplained; 

1.  The  measured  gain  is  larger  than  predicted  bv  Shields' 
theoretical  calculation. 

2.  In  10%  H?  -  90%  the  values  of  AT^  calculated  from  the  time 

dependence  of  the  AT  data  and  known  value  of  x  ,  x  and  t  are 

1  tc  vc  v 

unreasonable  when  compared  with  values  estimated  from  the  amount  of 
discharge  energy  going  into  vibration. 


CHAPTER  V 


CONCLUSIONS 

An  investigation  to  observe  sound  amolif ication  in  a  vibrationallv 
excited  gas  has  been  described.  Bv  means  of  an  electrical  discharge, 
nitrogen  gas  was  vibrationallv  excited  to  a  vibrational  temperature  of 
the  order  of  a  few  thousand  degrees  Kelvin.  The  discharge  was  also  able 
to  generate  a  sound  nulse  that  nropagated  through  the  excited  gas  after 
the  discharge  had  ceased.  Rv  discrete  Fourier  analysis  of  the  sound 
pulse  upon  each  reflection  in  the  discharge  tube  it  was  possible  to 
measure  the  absorption  and  velocitv  for  different  frequencies  as  a  func¬ 
tion  of  time.  Bv  comparing  the  measured  velocitv  with  the  numerically 
calculated  velocity,  the  variation  in  translational  temperature  of  the 
gas  with  time  could  be  determined.  In  some  cases  (50%  He  -  50%  at  20 
Torr,  50%  He  -  50%  at  40  torr  and  10%  -  90%  N^)  measureble  gain  was 

observed.  In  two  cases  (50%  He  -  50%  at  40  torr  and  10%  H^  -  90% 

at  40  Torr)  it  was  possible  to  measure  the  variation  of  gain  with  time. 

The  measured  decay  rate  for  the  difference  between  translational 

temperature  and  wall  temprature  (AT)  agreed  in  all  gas  mixtures  and 

pressures  with  values  calculated  from  independent Iv  determined  values  of 

the  vibrational  relaxation  time,  xv*  and  the  thermal  diffusion  time,  x  , 

assuming  that  the  thermal  diffusion  time  for  conducting  trans1 ational 

energy  to  the  wall  was  the  same  as  that  for  conducting  vibrational  energv 

to  the  wall  (x  =  x  )• 
tc  vc 
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Kxperiments  on  pure  nitrogen  confirmed  Nighan's  theoretical  calcula¬ 
tion^  from  measured  electron  cross  section  that  practiallv  all  the 
discharge  energy  goes  into  vibration  in  N  ,  i.e.,  the  discharge  was  found 
to  raise  the  translational  temperature  in  Dure  only  a  few  degrees.  In 
the  gas  mixtures,  ;he  discharge  was  found  to  raise  the  translational 
temperature  from  10  to  40°C  depending  upon  the  energy  per  mole  deposited. 

In  those  cases  where  the  gain  was  measured  it  was  about  an  order  of 
magnitude  greater  than  expected  from  Shields'  calculations.  These  calcu¬ 
lations  were  based  on  the  assumption  that  the  temperature  dependence  of 

the  transition  rate  was  given  by  -^-ln( x  *)  =  0.01.  Bv  increasing  this 

dT  v 

derivative  to  0.05  and  assumming  70%  of  the  discharge  energy  was  in 

vibration  at  t  ;  i  ,  Shields'  calculations  would  give  the  approximate 

tc 

magnitude  of  gain  in  both  the  gases  where  this  time  dependence  was  mea¬ 


sured  (50%  He  -  50%  M2  and  10%  H?  -  90%  at  40  Torr). 

In  the  10%  H„  -  90%  N_  gas  at  40  Torr,  t  /t  was  small  enough  so 
2  2  v  tc 

that  it  was  possible  to  imply  the  AT  vs  time  curve  from  AT  versus  time 

v 

curve.  However  AT  values  thus  obtained  indicate  that  only  about  8%  of 


v 

the  energy  stored  in  the  capacitor  was  coupled  into  the  vibrational 
energy  of  the  gas. 

The  time  averaged  gain  in  dB/cm  over  a  freauencv  range  from  1.9  to 
5.5  kHz  was  found,  within  the  accuracy  of  measurement,  to  be  independent 
of  frequency  as  predicted  bv  Shields'  calculations. 

Measurements  were  also  tried  in  a  20%  H?  -  80%  N?  mixture  at  a 
pressure  of  40  Torr.  The  signal  pattern  was  unrepeatable.  It  was 


observed  chat  even  though  the  first  reflections  of  the  sound  Dulse  for 

different  sparks  coincided,  the  following  reflections  were  randomlv 

shifted  along  the  time  axis  (on  the  oscilloscooe  screen).  This  meant 

that  the  velocity  (and  thus  the  gas  temperature)  was  not  varying  with 

time  in  the  same  way  from  soark  to  spark.  This  was  attributed  to  the 

small  value  of  the  ratio  x  lx  •  Evidently  the  gas  was  being  heated 

v  tc 

unevenly  by  the  vibrational  energy  relaxing  into  translation  faster  than 
it  could  be  conducted  to  the  wall.  Due  to  these  difficulties  the  experi¬ 
mental  data  for  20%  H^  -  80%  ^  mixture  was  not  reported  in  this  studv. 
However,  qual i tat i vel v ,  the  gain  was  observed  in  this  mixture  to  have  a 
slightly  larger  value  than  the  gain  observed  in  10%  -  80%  mixture. 

The  results  reported  here  suggest  two  definitive  experiments.  (1) 
Make  measurements  on  mixture  of  and  H^O.  The  H^O  will  shorten  the 
relaxation  time  but  should  produce  less  gain  since  the  vibrational  relax¬ 
ation  rate  is  nearly  temperature  independent.  (2)  A  second  experiment 
suggested  is  on  mixtures  of  CO  and  CO  like  should  be  left  highlv 

vibrationally  excited  by  the  discharge.  In  this  case  the  vibrational 
temperature  of  the  CO  can  be  monitored  by  the  infrared  emission  from 
vibrational  transition. 
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APPENDIX  A 


ComDuter  program  used  to  fourier  analyze  the  sound  nulse  and 
calculate  the  absorotion  and  velocitv. 


10  OPTION  BASE  1 

20  SHORT  V(256)  ,R ( 2 56 > ,1(256) 

30  SHORT  H(30,7),A(30),C(30,7),T(30),H1(30,7),T1(30,7) 

40  SHORT  S1(7),S2(7),T0(30,7),S3(7),S4(7),A1(30),C1(30,7) 

50  DIM  A$ [ 8200 ] ,BS[280] 

60  DISP  "GAS  SAMPLE" 

70  INPUT  GS 

80  DISP  "INPUT  P ,L ,T9" 

90  INPUT  P9.L9.T9 
100  PRINT  "GAS  SAMPLE'" ;GS 
110  PRINT  "PRESSURE'" ; P9 ; "MM" 

120  PRINT  "L  TUBE=";L9; "CM" 

130  PRINT  "R. TEMP. =" ; T9 ; "DEG. C" 

140  DISP  OF  POINTS  N,  //  OF  REFLECTION  Ml,  FREQ.  COMP  M2,  #  n F  SPARK 
M3" 

150  INPUT  N, Ml, M2, M3 

160  PRINT  "#  OF  SAMPLE  POINTS'" ; N 

170  T=N*. 00001 

180  FOR  M=1  TO  Ml 

190  FOR  1=1  TO  M2 

200  HI ( I ,M)=0  0  T1 ( I ,M)=0  3  C1(I,M)=0 

210  S3(M)=0  3  S4(M)=0  0  A1(I)=0 

220  NEXT  I  (3  NEXT  M 

230  FOR  K9=l  TO  M3 

240  PRINT  (3  PRINT  "SPARK  F\K9 

250  FOR  M=  1  TO  Ml 

260  DISP  "INPUT  S1(";M;") ,S2(";M;")" 

270  INPUT  SI (M) ,S2(M) 

280  S  3 ( M) =S  3 (M) +S 1 ( M ) 

290  S4(M)=S4(M)+S2(M) 

300  PRINT  "REFL  F';M; SI (M) ; S2(M) 

310  NEXT  M 

320  FOR  M= 1  TO  Ml 

330  GOSUB  1480 

340  FOR  I=S2(M)-S l(M)+2  TO  N 

350  V(I)=0 

360  NEXT  t 

370  FOR  1=1  TO  N/2 
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380  IF  2*1— 1=1  THEN  440 
390  IF  2*1=S2(M)-S1(M)+1  THEN  440 
400  IF  2*1=S2(M)-S1(M)+1  THEM  460 
410  R( I )=V ( 1*2-1 ) 

420  I(I)=V(I*2) 

430  GOTO  470 

440  R(I)=.5*V(2*I-1) 

450  GOTO  420 

460  I(I)=.5*V(2*I)  n  R( I ) =V (2*1—1 ) 

470  NEXT  I 
480  P=1 

490  FOR  1=1  TO  10 
500  P=P*2 

510  IF  P=N  THEN  530 

520  NEXT  I 

530  Pl=I 

540  N2=N/2 

550  J=0 

560  GOSUB  1100 

570  FOR  1=1  TO  M2 

530  TO ( I ,M)=ATN2 ( 1(1) ,R(I))-2*PI*(I-1)*S1(M)/N 
590  H(I ,M)=SQR( I(I)“2+R( I ) “2) 

600  T1(I,M)=T1(I ,M)+T0( I ,M) 

610  H1(I,M)=H1(I,M)+H(I,M) 

620  NEXT  I 

630  BEEP  (?  NEXT  M 

640  GOSUB  1670 

650  FOR  1=1  TO  M2 

660  A1(I)=A1(I)+A(I)  0  NEXT  I 

670  GOSUB  1850 

680  FOR  M= 1  TO  Ml-1 

690  FOR  1=1  TO  M2 

700  C1(I,M)=01(I,M)+C(I,M) 

710  NEXT  I  (?  NEXT  M 

720  IF  K9=M3  THEN  750 

730  BEEP  0  BEEP  <3  DISP  "INPUT  SPARK#" ;K9+1  0  PAUSE 

740  NEXT  K9 

750  BEEP  0  BEEP  (3  BEEP  0  piSP  "HONE!!!" 

760  FOR  1=1  TO  M2 
770  a1(I)=A1(I)/M3 

780  PRINT  I— 1 ; A1 (I)  (3  NEXT  I 
790  FOR  M= 1  TO  Ml-1 

800  PRINT  (3  PRINT  "VEL.  OF  REFL.  #"  ;M ; 

810  FOR  1=1  TO  M2 
820  C1(I,M)=C1(I ,M)/M3 

830  PRINT  1-1 ; C 1 ( I ,M) 

840  NEXT  I  13  NEXT  M 
850  FOR  M= 1  TO  Ml 
860  S3=S3(M)/M3 
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870  S4=S4(M)/M3 

880  PRINT  3  PRINT  "REFL. #" ;M; "STARTS  AT";S3;"ENDS  AT";S4 

890  FOR  1=1  TO  M2 

900  H(I.M)=H1(I,M)/M3 

910  T0(I,M)=T1(I,M)/M3 

920  PRINT  USING  930;  1-1 ;H(I ,M) ;TO( I ,M) 

930  IMAGE  40 , 2X ,MZ , 3DE , 2X ,MZ , 3DE 
940  NEXT  I 
950  NEXT  M 

960  PRINT  3  PRINT  "SLOPES  OF  LOG  ( H )  VS,M" 

970  GOSUB  1670 
980  FOR  1=1  TO  M2 
990  PRINT  I— 1 ; A( I)  (3  NEXT  I 
1000  GOSUB  1850 
1010  FOR  M=1  TO  Ml  —  1 
1020  PRINT  "VEL.  OF  REFL. #" ; M 
1030  FOR  1=1  TO  M2 
1040  PRINT  1-1 ;C( I ,M) 

1050  NEXT  I  @  NEXT  M 
1060  BEEP  3  BEEP  0  END 
1070  REM 

1080  REM  SUBROUTINE  1070  TO  1410  FOURIER  ANALYZES  DATA  POINTS 
1090  REM 
1 100  RAD 
1110  K=0 

1120  FOR  J=1  TO  N2-1 
1130  1=2 

1 140  IF  K<N2/I  THEN  1160 

1150  K=K-N2/I  3  1  =  1+1  3  GOTO  1140 

1160  K=K+N2/I 

1170  IF  K<=J  THEN  1190 

1180  A=R(  J+l  )  (3  R(  J+l  )=R(K+1 )  <3  R  (K+l  )  =A  0  p  =  I(.T+l)  (3  I  ( J+l  )  =  I(K+1 )  3 
I(K+1 )=A 
1190  NEXT  J 
1200  G=. 5  3  P2=l 

1210  FOR  1=1  TO  PI=1 

1220  G=G+G  (3  C=  1  3  E=0  3  0=SOR(  (  1=P2) /2) 

1230  P2=( l-2*( 1-1 ))*S0R( l+P2)/2) 

1240  FOR  R=1  TO  G 

1250  FOR  J=R  TO  N2  STEP  G+G 

1260  K=J+G  3  A=0*R(!<)+E*I(fO  3  R=F.*R(  K)-C*I  ( K)  (3  R(K)=R(J)-A  3 

I(K)  =  I(J)+B  3  R(J)=R(J)+A  (3  I(J)  =  I(J)-B 
1270  NEXT  J 

1280  A=F*P2+C*Q  3  C=C*P2-E*0  3  E=A 
1290  NEXT  R 
1300  NEXT  I 

1310  A=PT/N2  3  P2=COS  (A)  3  0=SIN(  A)  3  A=R(1)  3  R(n=A+l(l)  3  I(1)=A-I(1) 
1320  R( 1 )=R( 1 ) / 2  3  I(l)=l(I)/2 
1330  C=1  3  E  =0 
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1340  FOR  J=2  TO  N2/2 

1350  A=E*P2+C*0  *  C=C*P2-E*Q  0  E=A  «  K=N2-J+2  <?  A=R(J)+R(K) 

1360  B=(I(J)+I(  K)  )*C-(R(  J)-R(K)  )*E  (3  U=I(J)-I(K) 

1370  V=( I ( J )+I ( K) )*E+(R(J)-R(K))*C 

1380  R(J)=(A+3)/2  3  I(J)=(U-V)/2  <?  R(K)=( 4-R) /2  3  I (K ) =-( (U+V ) /2 ) 

1390  NEXT  J 

1400  I(  N2/2+1 ) =-I (N2/2+1 ) 

1410  FOR  J=1  TO  N2 

1420  R(J)=R(J)/N2  3  I(J)=I(J)/N2 

1430  NEXT  J 

1440  RETURN 

1450  REM 

1460  REM  SUBROUTINE  1420  TO  1560  TRANSFERS  DATA  POINTS  FROM  NICQLET  TO 
COMPUTER 
1470  REM 

1480  OUTPUT  715;  "NO" 

1490  IOBUFFER  BS 

1500  TRANSFER  714  TO  BS  FHS ;  COUNT  230 

1510  CLEAR  7 

1520  VO=VAL(BS  f  5 ,8]  ) 

1530  V1=VAL(BS[ 14,20]  ) 

1540  OUTPUT  715; "D3D2 
1550  IOBUFFER  AS 

1560  TRANSFER  714  TO  AS  FHS;  COUNT  8192 
1570  CLEAR  7 
1580  J=0 

1590  FOR  I=S 1 (M)*2+l  TO  S2(M)*2+1  STEP  2 
1600  J=J+1 

1610  V(  J)=V1*(256*NUM(A$[I] )+NUM( AS[ T  +  l 1 )-6 5536* (NUM( AS [ I ] )>  =  128) ) 

1620  NEXT  I 
1630  RETURN 
1640  REM 

1650  REM  SUBROUTINE  1580  TO  1720  CALCULATES  THE  SLOPE  OF  LOG ( AMP L. )  VS. 

THE  #  OF  REFLECTION 
1660  REM 

1670  FOP  1=1  TO  M2 
1680  X0=0  3  X1=0  3  Y0=0  3  S=0 

1690  FOR  M= 1  TO  Ml 
1700  Y=LOG(H( I ,M) ) 

1710  X=M 

1720  X2=X~2  3  S 1 =X*Y 

1730  X0=X0+X  3  XI =X 1+X2 

1740  Y0=Y0+Y  3  S=S+S 1 

1750  NEXT  M 
1760  D=X1*M1-X0~2 

1770  ni=s*Ml-XO*YO 
1780  D2=X l*YO-S*XO 
1790  A(I)=D1/D 

1800  NEXT  I 


1810  RETURN 
1820  REM 

1830  REM  SUBROIJNTINE  1730  TO  1800  CALCULTES  VELOCITY 
1840  REM 

1850  EOR  M= 1  TO  Ml  —  1 
1860  FOR  1=1  TO  M2 
1870  J=I-1 

1880  T ( I )=T0 ( I ,M+1 )-T0 ( I , M) 

1880  C(I,M)=-(2*PI*J*L9/T/T( I)) 

1900  NEXT  I 
1910  NEXT  M 
1920  RETURN 
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APPENDIX  R 

TABLE  24 

Physical  Properties  of  gases  study  at  30C °K  used  in  calculations. 


Thermal 


Cas 


Molecular 

'.•'eight 


Conduct ivi tv 
(cal  cm  'deg  ' sec  ' ) 


Viscositv 

,  -1  -K 

(gr  c"1  sec  ) 


\'9  28.013  2.5  602.6x10  7  1.7848x10 


25He-75N?  22.011  2.25  860.2x10  7  1.8542xin_A 

50He-50N  16.008  2. on  U70xin"7  1.0370v]o_i 
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